Dynamics of a Mathematical Hematopoietic Stem-Cell Population Model

We explore the bifurcations and dynamics of a scalar differential equation with a single constant delay which models the population of human hematopoietic stem cells in the bone marrow. One parameter continuation reveals that with a delay of just a few days, stable periodic dynamics can be generated of all periods from about one week up to one decade! The long period orbits seem to be generated by several mechanisms, one of which is a canard explosion, for which approximate the dynamics near the slow manifold. Two parameter continuation reveals parameter regions with even more exotic dynamics including quasi-periodic and phase-locked tori, and chaotic solutions. The panoply of dynamics that we find in the model demonstrates that instability in the stem cell dynamics can be sufficient to generate the rich behaviour seen in dynamic hematological diseases.